PASJ: Publ. Astron. Soc. Japan , 1-??, 
© 2011. Astronomical Society of Japan. 
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Abstract 

Particle acceleration in the inner ~ 200 pc of the Galaxy is discussed, where diffuse TeV 7-rays have 
been detected by the High Energy Stereoscopic System (HESS) observation. The diffuse 7-ray emission 
has a strong correlation with molecular clouds with large velocity dispersion, indicating the presence of 
turbulence. It is argued that the turbulence may contribute to the acceleration of cosmic rays via stochastic 
acceleration. The stochastic acceleration may energize cosmic-ray protons up to ~ 100 TeV and electrons 
to ~ 1 TeV in a relatively tenuous medium. The difference in the efficiency between protons and electrons 
supports the hadronic scenario of the diffuse TeV 7-ray emission. 
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1. INTRODUCTION 

The inner few hundreds parscc of the Galaxy or sim- 
ply the Galactic center (GC) is an outstanding source of 
high energy radiation including 7-ray, X-ray, and non- 
thermal radio emission (e.g., Aharonian et al., 2006; 
Mayer-Hasselwander et al., 1998; Koyama et al., 1996, 
2007; Muno et al., 2003; Law et al, 2008b,a). In addi- 
tion, there exists a strong concentration of dense molec- 
ular gas in the GC. This molecular feature, the Central 
Molecular Zone (CMZ), consists of three major clouds, 
Sgr A, Sgr B2, and Sgr C (e.g., Fukui et al, 1977; Scoville 
et al., 1975; Liszt & Spikcr, 1995), and has a total mass of 
3 x 10 7 M© (Dahmen et al., 1998), being characterized by 
an unusually large velocity dispersion of 15km/s (1<t), and 
high temperature, 30— 300 K (Huettemeister et al., 1993; 
Martin et al., 2004; Oka et al., 2005). The CMZ appears 
to consist of numerous dense clumps with a volume filling 
factor of ~ 0.2-0.5 (Sawada et al., 2001). 

Aharonian et al. (2006) reported the discovery of dif- 
fuse TeV 7-ray emission by the High Energy Stereoscopic 
System (HESS) from the GC. They found a correlation 
between the 7-ray emission and molecular clouds (MCs), 
which strongly indicates that the 7-rays arc radiated by 
cosmic-ray (CR) protons via decay of neutral pions pro- 
duced by nuclear collisions with ambient atoms. An im- 
portant suggestion by Wommer et al. (2008) is that the 
spatial profile of diffuse 7-rays cannot be explained by su- 
perposition of point sources. They argued that it is likely 
that CRs are accelerated in the whole inter-cloud medium 
(ICM), otherwise the observed 7-ray emission would peak 
at the injection points because of relatively fast energy 
loss of CR protons. 

Crocker et al. (2011a) has recently developed a detailed 
steady-state model for the CR environment in the GC. 
They assumed that there is a fast wind (> lOOkm/s ) 
blowing out from the GC, that transports most of the 



CR protons above and below the disk. Such a scenario 
seems natural in view of recently found Fermi bubbles 
(Su et al., 2010; Dobler et al., 2010), while so far there 
is no direct evidence for the wind. In any case, they 
found a reasonable fit to the broadband spectrum with 
gas densities around n ~ 10cm~ 3 , magnetic field strength 
around B ~ 100 /iG, wind speeds around V wm d ~ a few 
x lOOkm/s. The required energy input rate ~ 10 39 ergs _1 
is consistent with the conventional idea that supernova 
remnant (SNR) shocks convert the explosion energy into 
CR with a conversion efficiency of ~ 10% for a supernova 
rate of ~ 1/2500 yr. A possible problem in their mod- 
eling is, however, a pure power-law injection spectrum. 
Recent observations of SNRs in 7-ray band indicates that 
there is a break in the primary CR spectrum accelerated 
by a SNR. The break energy seems to depend on the age 
of SNRs, but can appear in GeV range for middle age 
(~ 10000 yr), and in TeV range for young SNRs (~ 1000 
yr) (e.g., Abdo et al., 2010; Aharonian et al, 2007) . It 
is argued that the spectral break may be produced as a 
result of the interaction between SNRs and MCs (Malkov 
et al., 2005). Therefore, it is natural to expect that a 
break will also appear in the CR spectrum injected by a 
supernova exploded in/around the CMZ. The 7-ray spec- 
trum obtained by the HESS observation, however, does 
not show any significant spectral breaks up to ~ 10 TeV, 
indicating that the primary CR spectrum extends above 
~ 100 TeV. If this is the case, we need a reacceleration 
mechanism of pre-existing CRs above the spectral break. 

It is important to note that the unusually large velocity 
dispersion measured in the GC indicates the presence of 
strong turbulence. Given the fact that the peak positions 
of velocity dispersion coincident with the TeV 7-ray peaks 
(Aharonian et al., 2006), we suggest an alternative particle 
acceleration mechanism due to the turbulence. Namely, 
we discuss a possible contribution of stochastic accelera- 
tion (also known as the second order Fermi acceleration) 
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by Alfvenic turbulence in the GC environment. We find 
that the stochastic acceleration process with parameters 
relevant to the GC can energize CR protons up to ~ 100 
TcV in ~ 10 5 yr in a relatively tenuous medium, while the 
maximum energy of electrons is limited to ~ 1 TeV due 
to strong synchrotron cooling, which favors the hadronic 
scenario of the diffuse TeV 7-ray emission. 

2. STOCHASTIC ACCELERATION MODEL 

The original idea of stochastic acceleration proposed by 
Fermi (1949) is to invoke random motions of magnetic 
clouds with typical velocities of ~ 10 km/s. Charged par- 
ticles colliding with the clouds gain energy due to the dif- 
ference between head-on and head-tail collision probabil- 
ities. It is well-known, however, that the acceleration is 
far too slow if we take a mean free path (or an inter-cloud 
distance) of ~ 1 pc. This picture is later elaborated by 
considering pitch-angle and momentum-space diffusion by 
Alfvenic turbulence. The Alfven wave is an incompress- 
ible, weakly damped normal mode of magnetized plasmas. 
Because of its nearly dissipationless feature, Alfven waves 
are believed to be ubiquitous. The essential difference 
from the Fermi's original idea is that the scattering mean 
free path is now determined by the power of the turbu- 
lence spectrum at the resonant wavenumber, but not by 
inter-could distances. Note that the stochastic accelera- 
tion is invoked in many other astrophysical environments 
as well, including solar corona, black-hole accretion disks, 
lobes of giant radio galaxies, and 7-ray bursts. 

In this study, we consider only spatially integrated CR 
spectra in which detailed geometry is assumed to be unim- 
portant. Therefore, wc arc essentially considering the CR 
distribution in a box. The escape of CRs from the box 
either due to diffusion or advection is thus treated by 
an energy-dependent escape term with a time scale of 
t osc , i.e., it is independent of energy for the advection- 
dominatcd escape, but may depend on energy for the 
diffusion-dominated escape. The CR diffusion-convection 
equation then reduced to (c.f., Petrosian & Liu, 2004) 
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where E is the kinetic energy, N(E,i) =Anp 2 J /(£,x,p)dx 
with f(t,x,p) being the distribution function, Q(E) is the 
total injection flux, r acc and ti oss represent typical acceler- 
ation and energy loss time scales, respectively. Note that 
in the above equation we use the approximation that the 
particle energy is highly relativistic. The important pa- 
rameter for the stochastic particle acceleration is the spa- 
tial diffusion coefficient K. Assuming elastic scattering of 
CRs in the frame moving with Alfven waves, we may esti- 
mate the acceleration time scale as r acc = (jj^f) — 

V A 

where va represents the Alfven speed (e.g., Skilling, 1975). 
Although it is possible to calculate the diffusion coefficient 
k directly by using quasi-linear theory, this theory is based 
on a number of assumptions which may not necessarily be 
valid in general. In the present study, we instead use the 
Bohm scaling k= ^£r g coc E, where c is the speed of light, 



and r g is the particle gyroradius. This expression indicates 
that the mean free path of the particle can be written by 
£r s , meaning that £ represents the efficiency of scattering 
and is sometimes called the "gyrofactor" . Numerical val- 
ues of £ estimated from the solar modulation of low-energy 
(< GeV) CRs in the heliosphere are of order 10—100, but 
may become as small as ~ 1 (the Bohm limit) in the vicin- 
ity of strong shocks. Given the large velocity dispersion 
measured in the GC, we anticipate that there exist re- 
gions in the GC where Alfvenic turbulence is so strong 
that £ becomes the same order as in the solar wind. The 
typical acceleration time scale can then be written using 
parameters relevant to the GC as 



r acc ~2.1xl0yr( ^ 
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The maximum attainable energy through this process 
may be estimated by equating this to characteristic loss 
or escape time scales. We are mainly interested in 
CRs with energies much higher than 1 GeV. For pro- 
tons of such high energies, the dominant loss process is 
due to nuclear collisions with ambient atoms with an 
energy-independent characteristic time scale of ti OSSiPP ~ 
3.7 x 10 5 yr (n/100 cm -3 ) -1 . In the case of electrons, 
the ionization, bremsstrahlung, inverse Compton, and 
synchrotron losses are potentially important. In the 
GC environment, since the magnetic field strength is 
estimated to be > 50 fiG, and most likely > 100 /iG 
(Crocker ct al., 2010, 2011b), the dominant loss pro- 
cess is the synchrotron cooling with ri oss sync ~ 1.3 x 
10 6 yr(^/GeV)- 1 (S/100/xG)- 2 unless the ambient den- 
sity is too high n > 10 3 cm -3 in which bremsstrahlung with 
Tioss.brcms — 2-4 x 10 5 yr (n/100 cm" 3 )' 1 becomes impor- 
tant. In addition, if the suggested fast wind (> 100 km/s) 
exists (Crocker et al., 2011a, b; Crocker & Aharonian, 
2011), a corresponding escape time scale of r csc — 4.1 x 
10 5 yr(K, ind /100km/s)" 1 ( /l / 4 2pc) should also be taken 
into account, where h is the height of the 7-ray emitting 
region. Alternatively, large scale flows over ~ 100 km/s 
found in this region (Fukui et al., 2006; Molinari et al., 
2011) may also contribute to the escape. Since the nature 
of the flows does not change our argument, we shall con- 
sider only a global wind blowing out from the disk in the 
following discussion. 

Combining these time scales, the estimated maximum 
attainable energy by the stochastic acceleration as a 
function of ambient density is shown in Figure 1. In 
most cases, the maximum energy for protons is much 
larger than electrons because of strong synchrotron cool- 
ing which dominates at n ^ 10 3 cm -3 . For protons, the 
dominant loss changes from hadronic collisions at high 
density regions n > 100 cm -3 to the convective escape by 
the hypothetical global outflow in lower-density regions. 
Even in the presence of a fast wind of a few hundreds 
km/s that sweeps CRs out from the system in ~ 10 5 yr, 
it is possible to accelerate protons to ~ 100 TeV or even 
more when n ^ 10 cm -3 , that is lower than the volumetric 
average density n ~ 100 cm -3 in that region. It is im- 
portant to mention that the GC environment is highly 
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inhomogeneous, where MCs with much higher density 
10 3 cm -3 ) and tenuous ICM coexist. The modeling re- 
sult of Crocker et al. (2011a) actually favors the target gas 
density of n ~ 10cm~ 3 , suggesting that CRs accelerated 
in low-density regions do not penetrate deep into MCs. 
Therefore, if 7-rays are originating from such regions, the 
stochastic acceleration can provide a non-negligible con- 
tribution to the 7-ray flux. 

We now investigate the allowed parameter range of the 
gyrofactor £ in which stochastic acceleration alone is ca- 
pable of accelerating CR protons up to ~ 100 TcV; the 
constraint imposed by the HESS observation. This may 
have an important implication because it is generally diffi- 
cult to calculate the diffusion coefficient both theoretically 
and experimentally in a realistic, strongly turbulent envi- 
ronment. For this purpose, we may rewrite equation (2) 
by equating r acc with an unspecified loss or escape time 
scale r, 



case, the resulting steady-state spectrum above the injec- 
tion energy (Q(E) = 0) with negligible losses (ti oss = 00) 
becomes a power law. More specifically, by substituting 
these into equation (1) and assuming a power-law func- 
tion N(E) oc E~ p , it is easy to show that the power-law 

index p should satisfy p = — 1 + ^/1 + 5^. The ratio of two 
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Since nuclear collisions are not the dominant loss in low- 
density regions n < 10cm -3 , r is essentially determined by 
the wind velocity (Crocker et al., 2011a). A wind velocity 
of a few hundreds km/s indicates that £ ~ 10—100 is the 
necessary condition for B = 100— 300 /iG. These values are 
similar to those estimated in the heliosphere. We thus ar- 
gue that there exist regions in the GC where such strong 
diffusion and corresponding stochastic acceleration are ac- 
tually taking place. A possible example of such scenario 
is HESS J1745-303, which has been discussed recently by 
Hayakawa et al. (2011). 

The fact that the GC is far from uniform, and the CR 
production and radiation may occur in different environ- 
ment indicates the importance of understanding the in- 
teractions between these regions. In general, the accel- 
eration region may have a finite bulk velocity relative to 
MCs. It is therefore natural to consider collisions between 
the two distinct regions. Actually, we expect stronger 
hadronic 7-ray emission from such an interaction region 
than more tenuous (i.e., not interacting with MCs), but 
potentially more efficient acceleration regions. The con- 
vective transport time scale of CRs deep into a MC is 
estimated as ~ 6.5 x 10 5 yr(//10pc)(cr/15km/s)~ 1 , where 
I is the typical size of a MC, and a is the velocity disper- 
sion. If CRs were merely convected by the background 
flow, they would lose energies faster than this time scale 
because the hadronic-collision loss becomes increasingly 
significant during the convective transport. In reality, 
CRs can diffuse against the flow and thus remain in the 
acceleration region longer than expected above. In the 
following, we consider only the case of a one-dimensional 
interaction for simplicity. The diffusion length of CRs 
against a bulk flow u can be written as k/u. We may 
thus estimate the time scale of CR escape to a dense 
MC as r osc ~ k/u 2 , where the acceleration is essentially 
terminated. It is important to notice that r acc and T esc 
have the same energy dependence, i.e., both are propor- 
tional to oc E under the Bohm scaling assumption. In this 



timescales T acc /T csc ~ (u/va) 2 thus determines the power- 
law index. A measured spectral index of ~ 2.3 requires a 
fast colliding velocity of u ~ ?>.1va which is similar to the 
inferred wind velocity of a few hundreds km/s. It is pos- 
sible to interpret that the measured 7-rays are originating 
mainly from such strongly-interacting regions, although 
we think that the theory at the present stage may over- 
simplify the situation for quantitative comparisons. This 
is partly because we assume that the acceleration region is 
a homogeneous medium, while interactions with MCs in- 
evitably produce inhomogeneity, particularly gradients in 
density and magnetic field. Other reasons include steady- 
state, and one-dimensional assumptions, etc. Therefore, 
the present model must be taken as a qualitative model 
rather than a quantitative one. Nevertheless, we consider 
that the stochastic acceleration can in principle play a role 
for accelerating > TeV energy CR protons and electrons 
in the GC, particularly in the low-density ICM. 

3. SUMMARY AND DISCUSSION 

Aharonian et al. (2006) have proposed a time-dependent 
scenario in which the diffuse CRs in the GC are pro- 
duced by a putative past event occurred ~ 10 4 yr ago, and 
then diffuse over a ~ 100 pc scale. Recent picture of the 
GC is, however, such that particle acceleration to multi- 
TeV energies occurs not only in individual (i.e., point-like) 
sources but also in the whole ICM. This is because sim- 
ple spatial diffusion from a single source will produce a 
peak at the source position, which is not consistent with 
the observed spatial profile. Also, no evidence for diffusion 
hardening due to an energy-dependent diffusion coefficient 
at the propagation front is reported so far. 

An important indication from observations is that large 
velocity dispersion comparable to Alfven velocity exists 
throughout the GC, indicating the presence of turbulent 
flows. The total turbulent kinetic energy may be esti- 
mated as ~ 10 53 erg, while that of diffuse CR component 
is ~ 10 49 erg in 4—40 TeV and ~ 10 50 erg if the spectrum 
extends from 10 9 to 10 15 eV (Aharonian et al., 2006). 
Therefore, the observed turbulence has sufficient energy 
to produce the diffuse CRs whatever the mechanism is. 
In the present study, we argue that the stochastic accel- 
eration by Alfvenic turbulence can accelerate CRs up to 
^100 TeV for protons, and ~ 1 TeV for electrons within a 
reasonable parameter range. The difference in the maxi- 
mum energy between two components favors the hadronic 
scenario of the TeV 7-ray emission, consistent with the 
correlation found between MCs and 7-rays. Of particu- 
lar note is that the GC is highly inhomogeneous; dense, 
but clumpy MCs are embedded in much more tenuous, 
volume-filling regions where the stochastic acceleration is 
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more efficient. 

Crocker et al. (2011a) has also suggested that the 7-ray 
emitting region is less dense ~ 10cm~ 3 than the volumetric 
average density of ~ 100 cm -3 , which is basically consis- 
tent with our idea. Our theory is compatible with their 
assumption that the wind blowing from the GC trans- 
ports most of the CR protons that illuminate the Fermi 
bubbles. However, the lower limit of the wind velocity de- 
termined to explain radio synchrotron emission from the 
GC lobe does not apply in our case. This is because the 
stochastic acceleration can overcome synchrotron cooling 
in the energy range of radio-emitting electrons. If an as- 
sumed empirical relation between the total infrared and 
radio-continuum (and also 7-ray) luminosities applies to 
the GC, a fair lower limit would be ~ lOOkm/s to sweep 
out CR protons before they are convected into the cores 
of MCs (Crocker et al., 2011b). In any case, it is likely 
that CRs accelerated in the ICM emit hadronic 7-rays by 
interacting with MCs. With this in mind, we investigate 
collisional interactions between the acceleration regions 
and MCs that may produce a power-law type spectrum. 
More detailed studies on this subject are still needed. 

We also comment on a similar work done by Melia & 
Fatuzzo (2011). While their idea itself is essentially the 
same as the stochastic acceleration model, their equation 
for the energy gain is too optimistic as also noted by them- 
selves. Their model at least differs from the standard 
stochastic acceleration model, which must be considered 
with great care. 

Finally, we would like to mention that the stochastic 
acceleration is of course not the only candidate mechanism 
of CR acceleration in the GC. Magnetic reconnection and 
shocks (not associated with SNRs) may also contribute 
acceleration of CRs in the turbulent GC (e.g., Lesch & 
Reich, 1992), which should be elaborated in the future. 

This work is in part supported by the Global COE 
Program of Nagoya University (QFPU) and KAKENHI 
22740118 (T.A.) from JSPS and MEXT of Japan. 
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